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Abstract 

A variational approach is developed for bound state calculations in three- and four-electron 
atomic systems. This approach can be applied to determine, in principle, an arbitrary bound 
state in three- and four-electron ions and atoms. Our variational wave functions are constructed 
from four- and five-body gaussoids which depend upon the six (ri2, ri3, ri4, r23, r24, 7*34) and ten 
(ri2, ri3, ru, r%s, r^, r24, r25, r^, r^ and r^) relative coordinates, respectively. The approach al- 
lows one to operate with the different number of electron spin functions. In particular, the trial wave 
functions for the ^-states in four-electron atomic systems include the two independent spin func- 
tions xi = ot(3al3+f3a(3oL — fioLafi — aPfia and \2 = 2aa(3(3+2(3(3aa — (3aa(3 — af3f3a — (3a(3a — a(3af3. 
We also discuss the construction of variational wave functions for the excited 2 3 S"-states in four- 
electron atomic systems. 
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I. INTRODUCTION 

In this study we consider the electronic structure of three- and four-electron atoms and 
ions. In particular, we undertake variational computations of the singlet ^-states and triplet 
3 S— states in various four-electron atoms and ions. Below, by an atomic system we mean a 
system which contains a number of electrons and one heavy nucleus. In the non-relativistic 
approximation used in this study the Hamiltonian of an arbitrary (A — l)-electron atomic 
system takes the form (see, e.g., 



where A is the total number of bodies in the atomic system. In a three-electron atomic system 
one has A = 4 in Eq.flT]). In this case the subscripts 1, 2, 3 designate three electrons, while 
the subscript 4 denotes the positively charged atomic nucleus. For four-electron systems 
we have A = 5 in Eq.flJJ and the subscripts 1, 2, 3, 4 designate four electrons, while the 
subscript 5 denotes the positively charged nucleus. Note that the Hamiltonian, Eq.flJJ), and 
all equations that follow are written in atomic units, where h = l,m e = l,e = 1. Also, 
in Eq.flTJ Vj = ^-) is the gradient operator of the i— th particle {% = 1, 2, . . . , A). 

The notation stands for the (ij)— relative distance/coordinate between i— th and j— th 
particles, i.e. =\ Ti — tj |= rji, where are the Cartesian coordinates of the z— th 
particle. Also, in Eq.flT]) the notation M stands for the mass of the central (heavy) nucleus, 
i.e. M > 1. 

The main goal of this work is to discuss some important details of bound state calculations 
of three- and four-electron atoms and ions. In particular, special attention will be given to 
the correct symmetrization of the four-electron trial wave functions which include more 
than one independent (electron) spin function, something not considered elsewhere in the 
modern literature. Formally, our main goal is to determine highly accurate solutions of 
the Schrodinger equation Hty = E^fr where H is the Hamiltonian, Eq.flT]), while E is the 
eigenvalue (= total energy) of the considered bound state and \I/ is the corresponding wave 
function. It is clear that the permutation symmetry of the total wave function \1/ must 
be different for the singlet and triplet bound states in four-electron system. The explicit 
construction of trial wave functions with the correct permutation symmetry between all 
three or all four electrons is the principal part of any accurate variational calculation of such 
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atomic systems. 

In general, bound state computations of three-, four- and many-electron atomic systems 
with the use of a number of different spin functions (which, however, correspond to the 
same value of the total (electron) spin S e and its z— projection (S e ) z ) are significantly more 
complicated than analogous calculations for two-electron helium-like atoms and ions. On the 
other hand, there is an obvious similarity in calculations of three- and four-electron atomic 
systems. Indeed, for the ground (doublet) l 2 S'-states in three-electron atoms and ions one 
needs to use the two different spin functions. The same number of spin functions is needed 
for accurate computations of the ^-states in four-electron atomic systems. 

Note that in any four-electron atom and/or ion all bound states are separated into two 
series of states: singlet states with the total electron spin S e equals zero, i.e. S e = 0; and 
triplet states with the total electron spin S e equals unity, S e = 1. In all previous works 
only the ground singlet l 1 ^— state was considered. Accurate computations of the triplet 
states in four-electron atomic systems have been performed for a very few ions/atoms 
[sj]. Moreover, in almost all previous computations of the singlet state in four-electron 
(Be-like) atoms and ions only one spin function xi — ot/3a/3 + (3a(3a — f3aa(3 — a/3(3a was 
used. Here and everywhere below in this study a and j3 are the spin-up and spin-down 
single electron functions, i.e. a z a = \a and a z (3 = — ~/3. The second independent spin 
function \2 = 2aa/3/3 + 2j3f3aa — f3aaj3 — aj3f3a — f3a(3a — af3af3 has been ignored in almost 
all modern accurate computations of the singlet states in four-electron atomic systems. 
Bearing this in mind we wanted to develop the method which can be used to perform bound 
state computations for the 1 5'— singlet and 3 S— triplet bound states in arbitrary four-electron 
atoms and ions. Our method is not restricted with respect to the number of spin functions 
included. It works equally well in cases, when one, two, three and even more independent 
spin functions are used. 

II. VARIATIONAL WAVE FUNCTIONS 

A central feature of any variational method is the construction of trial wave functions \& 
with the correct permutation symmetry. In general, such a trial wave function must include 
all electron and nuclear coordinates. Accurate wave functions explicitly depend upon all 
scalar interparticle coordinates =| r— r,- | in the atomic system. The use of a large number 
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of relative coordinates complicates the explicit symmetrization of the trial wave functions. 
Another complication follows from the presence of different (independent) terms in the spin 
part of the total wave function. For instance, let us discuss the case of the singlet state 
in four-electron atomic systems. To compute this state in this study we shall use the two 
following independent spin functions written in the form \i — otj3aj3 + (3a(3a — [3aa(3 — a(3f3a 
and X2 = 2aa(3{3 + 2f3(3aa — j3aa(5 — af3(3a — f3aj3a — a{3a(3. These two functions obey the 
following relations 

S 2 Xk = , (S z ) Xk = , (Xi\Xk) = D k 5 ik (2) 

where k = 1,2; D k are the normalization factors of the spin functions and S = s 4 + s 2 + s 3 + s 4 
is the total electron spin of the four-electron system. 

The total wave function of the 1 S(L = 0)-state of the four-electron beryllium-like 

n n 

atom/ion is written in the form (see, e.g., [4|, [5[ and references therein) 

^l=o = ^l=o(A; {ry}) • X i + <f>L=o(B; {r %] }) ■ X 2 (3) 

where if)L =0 (A; {r^}) and </)l=o(B; {r^}) are the two independent radial parts (= spatial 
parts) of the total wave function. For the wave function, Eq.()3]), one finds S 2 \& = and 
S z ^> = 0, where S = s x + s 2 + s 3 + s 4 is the total electron spin of the four-electron system. 
The radial parts if)L =0 (A; {rij}) and (J)l=o(B; {r^}) of the total wave function Eq.([3]) in this 
work are represented in the form {(| 

N A 

i> L=Q (A; {r y }) =VxY,C k - «p(- E 4H) ( 4 ) 

k=l ij 

N B 

<Pl=o(B; {nj}) =V 2 J2C k - ex P (- £ /f Vj) (5) 

k=l ij 

where Na and Nb are the number of basis functions used, C k and C k are the linear parameters 
of the variational expansions, Eq.(j4]) and Eq.(j5]), while {r^} is the set of relative coordinates 
which are needed for complete description of five-body systems. The notations a[j and 

(k) 

plj designate the non-linear parameters associated with the relative coordinate in the 
k— th basis function in the i/j and 4> radial functions, respectively (see, Eqs.()l]) - §5§). For 
all beryllium-like ions and atoms considered in this study the notation {r^} stands for ten 
relative coordinates ri 2 , ri 3 , ri 4 , ri 5 , r 2 3, r 24 , r 25 , r 34 , r 35 and r 45 . The radial basis functions in 
Eqs.Q - ([5]) are called the five-body gaussoids of the ten relative coordinates. This name 



was used in [6] where these basis functions were invented for nuclear few-body systems. 



The main advantage of the radial function defined in Eqs.(jll) - Eqs. ([3]) follows from the 
fact that the formulas for all matrix elements do not depend explicitly upon the total number 
of particles in the system. In other words, these formulas are essentially the same for three-, 
four-, five- and many-body systems and are discussed in the next Section. 

The notations A and B in Eq.fll]) and Eq.([I)]) mean that there are two different sets of non- 
linear parameters. Note that each of the basis functions in Eq.Q and Eq.Q contains ten 
non- linear parameters which are optimized independently. The summation over (ij) = (ji) 
in the exponents of Eq.(j3]) and Eq.([5]) is taken over all possible different pairs of particles. 
In general, the radial basis functions are not orthogonal to each other. The projectors V\ 
and Vi produce the trial wave functions with the correct permutation symmetry between 
all electrons (see below). The symbol L in Eq.fllj) is used for the total angular momentum 
of the considered system. For the ground state of any Be-like system we always have L = 
and the total electron spin of such states equals zero. Furthermore, the parity of these states 
in the four-electron systems is even. With respect to this, these states are often designated 
as the 1 S e states, or l 1 ^ 6 states. 

For the triplet 3 S— states in four-electron atomic systems one also finds the two inde- 
pendent spin functions Xi — ct(3aa — f3aaa and X2 = 2aa(3a — (3aaa — af3aa. Note that 
there are, in fact, the two independent triplets of the spin functions, i.e. the total number 
of electron spin functions equals six (6 = 2 x 3). The functions Xi an d Xi mentioned above 
correspond to the values S e — 1 and (S e ) z = 1. By using the explicit form of these two spin 
functions one can construct the four remaining spin functions with S e = 1 and (S e ) z = 
(two functions) and S e = 1 and (S e ) z = —1 (two functions). In actual calculations of inter- 
nal atomic structure we can restrict ourselves to the use of the two spin functions Xi an d 
X2 only. All six spin functions are needed only in some special cases, e.g., if an external 
magnetic field is present. 

It is interesting to find that the total variational wave function \l/ of three-electron atoms 
and ions is also represented in the form of Eq.flBJ). In this case xi — a P a ~ fi aa an d 
X2 = 2aa/3 — j3aa — a(3a, while the two independent radial parts depend upon six relative 
coordinates r 12 , r 13 , r 2 3, r 14 , r 2 4, r 34 (here the indexes 1, 2, 3 mean three electrons, while 4 
stands for the nucleus). The radial basis functions can also be chosen in the form of Eq.()2]) 
and Eq.([5]), i.e. in the form of four-body gaussoids. In this case each of the radial basis 
functions contains six non-linear parameters. Note also that for each of the spin functions 
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we have S 2 Xi = \x% an d (S)zXi = (here 2 = 1, 2). The explicit constructions of the 

variational wave functions for three- and four-electron atomic systems is discussed below. 



III. MATRIX ELEMENTS 

Actual computation of the matrix elements with the four- and five-body gaussoids Eq.flU) 
and Eq.flHJ), is based on analytical formulas derived elsewhere (see, e.g., jf|, Q). It was 
mentioned already in |(| that the explicit expressions for all matrix elements needed for 
the solution of the Schrodinger equation depend upon the total number of particles A in 
the system as a numerical parameter. In other words, the expressions for matrix elements 
in three-electron systems coincide with the corresponding formulas for matrix elements ob- 
tained for four- and five- and many-electron systems. Moreover, in some few-body systems 
one of the electrons can be replaced by another particle, e.g., by \i~ muon, but such a re- 
placement does not change the explicit formulas for matrix elements. Below, the symbol A 
designates the total number of particles (i.e. bodies) in the system. The A— particle atomic 
system includes the (A — 1)— electron subsystem plus one heavy nucleus. The A— particle 
muonic atom/ion contains the (A — 2)— electron subsystem, one negatively charged muon 
fi~ and one heavy nucleus. 

The explicit formulas for all matrix elements needed in computations of the A-particle 
atomic systems can be presented in a very brief form with the use of the following compact 
notation 



(a \= (a« |= exp(- £ a% • rj) and | 0) =| ^) = exp(- £ % ■ r%) (6) 



where A is the total number of particles in the system. In this notation the symbols | (3) 
and | a) (or (a \ and (/3 |) stand for the radial basis functions. These notations are different 
from the symbols a and (3 used in other Sections of this study to designate the spin-up and 
spin- down functions. 

In the notation defined in Eq.Q the overlap matrix element (a \ j3) is written in the 
form 



.4 



A 



i>j=l 



i>j=l 



(a\P) = (a {k) | p®) 




(7) 



= 71 
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where D is the determinant of the (A — 1) x (A — 1)— matrix B with the matrix elements 

^ = XX+4), i?j = l,2,...,A-l (8) 

= -(«&• + #;), i?j = l,2,...,A-l 

In particular, the explicit expression for the (k,£) matrix element of the overlap matrix S 
in the case of A = 5 is the 4x4 matrix B with matrix elements 6^- defined in Eq.flSJ). 
Analytical and/or numerical computations of the determinant of this matrix and all its first 
order derivatives is straightforward. 

The formula for the appropriate matrix elements of the potential energy can be written 
in the form: 



A r+oo , I ) ■ \ 

5> | V(r tJ ) | P) = —(a \ 0)<£ V N^) • exp(-x 2 ) • x 2 rfx (9) 
(ij) V ij 

where D tj = j£- = while {ij) = (ji) = (12), (13), (23), (14), (24), (34) for the three- 
electron atom/ion (A = 4) and (ij) = (ji) = (12), (13), (14), (15), . . ., (35), (45) for the 
four-electron atom/ion (A = 5). The explicit expressions for various interparticle potentials 
often used in bound state calculations can be found in [6]. The integral in the last formula 
is computed analytically in many actual cases, including the case of Coulomb, Yukawa-type, 
exponential, oscillator, and other potentials. The matrix elements of the kinetic energy take 
the form (in atomic units) 

(« l T l^> = ^[E —(A* + D jk - D l3 )] (a | 0) (10) 

lD l ijk =l m k 

where rrii(i = 1,2,..., A) are the masses of the particles and i ^ j ^ k in Eq. ffTO]) . 
The explicit formulas for matrix elements of other operators written in the basis of many- 

n n 

dimensional gaussoids can be found elsewhere (see, e.g., |6(, [7|). 

If all formulas needed for matrix elements of the potential and kinetic energies are known, 
then the solution of the incident Schrodinger equation is reduced to the following generalized 
eigenvalue problem 

N 

Y,{ H ^-E- s a ,e)Cp = o (n) 
0=1 

for a = 1,...,N, where iV is the total number of basis functions used. Here H a p = 
T a ,/3 + V a $ — (oe | T | /?) + {a | V \ j3) is the Hamiltonian matrix, while T Q;( g = (a \ T \ (3) 



and V aj p — (a \ V \ j5) are the matrices of the kinetic and potential energies, respectively. 
The S a> p = (a \ (3) matrix in Eq. fllll) is the overlap matrix, Eq.([7]). For non-orthogonal 
basis sets the overlap matrix is a typical dense matrix, i.e. all elements of such a matrix 
differ, in general, from zero. Moreover, it can be shown that the overlap matrix (a \ j3) is a 
symmetric, positively defined matrix. This means that all eigenvalues of the overlap matrix 
are positive. 

IV. ANTISYMMETRIZATION OF THE TRIAL WAVE FUNCTIONS 

Let us consider the antisymmetrization of the trial wave functions and related anti- 
symmetrization of the corresponding matrix elements derived in the previous Section. As 
mentioned above the correct antisymmetrization is a central part of the construction of 
explicitly correlated, trial wave functions. In general, such a wave function depends upon 
all electron-nuclear and electron-electron coordinates. In two-electron atoms and ions the 
antisymmetrization of the total wave function is a trivial problem, since the wave function 
of the two-electron system is always represented as a product of a radial and two-electron 
spin functions. Moreover, only singlet and triplet spin functions are possible in any two- 
electron atom and/or ion. The singlet states have spin function xi — a P — flat, while triplet 
states have three spin functions X2 = aot -> X2 = a P + P a -> X2 = PP- F° r the singlet 
spin function one finds S 2 Xi = 0, S z Xi = 0, while for the triplet spin functions we have 
S 2 X2 = 1(1 + 1)X2 ] = 2 X2 } and S Z X2 = K iXi\ where k x = 1, k 2 = and k 3 = -1, 
respectively. Here and below S and S z are the total electron spin of the system and its 
z— projection. 

Note that the singlet Xi s pi n function is antisymmetric upon the electron variables. There- 
fore, its product with a symmetric radial function produces a function which is completely 
antisymmetric upon all electron variables. It is clear that such a function can be considered 
as a total wave function with the correct permutation symmetry between two electrons. For 
triplet states the corresponding radial function must be antisymmetric upon all electron 
(spatial) coordinates. 

In contrast with two-electron systems, the antisymmetrization of three-electron wave 
functions is a significantly more complex process, since for the same spin state one finds 
not one, but a number of different independent spin functions. This statement is true 
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even for the doublet states with S = | (i.e. 2S + 1 = 2) and S z = ±~ in any three- 
electron atomic system. In actual computations, such spin functions are usually chosen to 
be orthogonal to each other. For instance, any variational expansion written for the doublet 
^-states in a three-electron atomic system must include the two independent spin functions 
Xi = oi/3a — f3aa and xi — 2ao.fi — j3aa — a/3a. The total wave function for the ground 



doub 
e.g., 



ct 
9]. 



S(L = 0)— state of the three-electron atomic system is written in the form (see, 



id) 



\1/ L=0 = ifj L=0 (A; {r^ }){aj3 a — j3aa) + <pL = o(B; {rij})(2aaj3 — j3aa — af3a) (12) 

where if)L =Q (A; {?>,•}) and 0l =o (-E>; {rij}) are the two independent spatial parts (= radial 
parts) of the total wave function. The notations A and B mean that the two sets of non- 
linear parameters associated with ip and (ft are optimized independently. 

Note that for each of these two spin functions Xi an d X2 the two following equations are 
obeyed: 

3 1 

S 2 Xfc = S(S + l)xk = -Xk , (S)zXk = ~Xk (13) 

where k = 1, 2 and S = s x + s 2 + s 3 is the total electron spin of the three-electron system, 
while S 2 is its z— projection. The conditions Eq.ffTBl indicate clearly that the two spin 
functions Xi an d X2 are equally important in this variational method. Therefore, we cannot 
neglect any of these spin functions a priori. This means that in any of our calculations for 
three-electron atomic system we have to make appropriate use of the two different radial 
functions and two spin configurations Xi an d X2- The explicit construction of the trial wave 
functions for three-electron atomic systems with two independent spin functions is more 
complicated problem than in the case of one spin function and its solution is based on the 
method of projection operators is discussed below. 



A. Three-electron atomic systems. Doublet states. 

Suppose our trial wave function for a three-electron atomic system is written in the form 
of Eq. ( TT2l) . In real applications, however, only those trial functions are accepted which have 
the correct permutation symmetry between all identical particles, i.e. electrons. This means 
that the two terms in the right-hand side of Eq. (|T2|) must be completely antisymmetric upon 
spin and spatial coordinates of the three electrons, i.e. upon the indexes 1, 2 and 3 in our 
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notation. In other words, we must have A e ^ = —*ff, where \& is given by Eq.( ll2p and A e is 
the three-particle (= electron) antisymmetrizer [ll|, 12] 



X = -(e - Aa - As - A3 + A23 + A32) (14) 
o 

Here e is the identity permutation, while Aj is the permutation of the z-th and j-th particles. 
Analogously, the operator Ajfc is the permutation of the i-th, j-th and k-th particles. The 
same notations are used everywhere below in the text. 

By using the three-particle antisymmetrizer, Eq.( !T4l) . we can construct a trial wave func- 
tion of the correct permutation symmetry. In reality, we need the matrix elements with the 
correct permutation symmetry, rather than the wave function itself. Let us describe the 
approach which allows one to obtain properly symmetrized matrix elements. First, note 
that the expectation value of an arbitrary completely symmetric operator W is written in 
the form 

(AE^(A; {rij})xi I W I A e J2^Mf, {rij})Xj) ( 15 ) 
where Xi are the spin functions (i = 1, . . . , N s ), where the A e operator is defined in Eq.(fl4l). 
The spin functions Xi an d Xi are assumed to be orthogonal to each other, i.e. (xi \ Xj) — 
The notations ipi(Af, {r^}) stand for the corresponding radial functions which depend upon 
all relative coordinates {r^} and non-linear parameters A±. These radial functions can be 
arbitrary, i.e. they are not necessarily orthogonal to each other. The operator W is a 
differential operator written in the relative coordinates. It is assumed to be completely 
symmetric in respect to all inter-electron permutations. 

Now note that completely symmetric operator W commutes with the A e operator, 



Eq.( TH|) . Moreover, the operator A e is an orthogonal projector [13|, i.e. (A e ) 2 = A e and 
(Ae)* — A e , where the symbol B* means the operator conjugate to the operator B. If the 
operator W does not depend upon spin variables, then by using these properties of the A e 
operator one can reduce the formula, Eq.( ll5j) . to the following form 

EE(^; {*■«}) I w I {r«}))txi I X I x 3 ) (i6) 

i 3 

This expectation value can be re-written into another form with the use of the following 
matrix notation 

mW)] ij = (xi\A e \Xj)[wMj (17) 
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where [. . ]ij designates the (ij)— matrix element of the corresponding matrix. The matrix 
elements of this operator (Q(W)) computed on any basis set of spatial three-electron wave 
functions have the correct permutation symmetry between all identical particles (electrons). 
The dimension of the [©(W 7 ")]^- matrix equals to the number of spin functions used in calcu- 
lations. 

First, let us compute the matrix elements of the A e operator, where A e designates the 
complete three-particle antisymmetrizer, Eq. (I14[) . Based on Eq.( lT7|) this operator can be 
written in the form 0(e), where e is the unit operator. Note that operator A e , Eq. (ll4p . is 
written in the form 

A e = S a b c Pabc (18) 
abc 

where s a b c is an integer number, while P a b c are the interparticle permutations in the system of 
three identical particles (see Eq.f JT4|) ). The sum in Eq.f JT8|) is computed over all interparticle 
permutations possible in three-body systems, Eq.( TH|) . The (ij)— matrix element of the A e 
operator in our basis is 

[*^-e]y ^ ^ ^abc(Xi \ Pabc \ Xj) Pabc (1^) 
abc 

In the case of the (ij)— matrix element of the Q(W) operator one finds 

[9{W)]ij = E S *bc(Xi I Pabc | Xi)WP abc = W[J2 Sabc(Xi I Pabc \ Xj) Pabc] (20) 
abc abc 

Note that both the expectation value (x% \ P a bc I Xj) an d s a bc are the integers for all abc, 
while the operator P a b c is a projector which acts on the spatial coordinates of three electrons 
(a — > b — > c). In other words, the A e operator is represented as the finite sum of all 
spatial permutations P a b c with integer coefficients equal to the products of s a b c , Eq.( fT5]) . 
and the (xi \ Pabc \ Xj) expectation values. The computation of all expectation values 
(Xi I Pabc I Xj) — (Xi I PabcXj) can De considered as the integration over electron spin 
coordinates. 

Based on Eq.f l20p we can introduce the following operator V 

V = V^ SabciXi I Pabc | Xj)Pabc (21) 
abc 

where T> is the normalization constant. The numerical value of D in Eq. (!2TT) is determined 
by the idempotiency of the V operator, i.e. V 2 = V. The explicit use of this operator 
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substantially simplifies all following formulas. For instance, the (a, /?) — matrix element of 
any arbitrary completely symmetric operator W can be written in the form 

(Va\W\ VP) = (a | VWV \ /3) = (a \ WV \ (3) = (a \ W \ V/3) (22) 

where a and (3 are non-symmetric basis functions. This matrix element has the correct 
permutation symmetry between all electrons. Moreover, all expectation values constructed 
with the use of these matrix elements, Eq. fl22p . also have the correct permutation symmetry 
between all electrons. This is the main advantage of constructing the orthogonal spatial 
projector V in explicit form. 

In actual computations of the doublet 2 5-states in three-electron atomic systems after 
the integration over electron spin coordinates, one finds the four following spatial projectors 

= ^= (2e + 2 P 12 - As - As - A23 - A32) (23) 

V^ = l - (As - A3 + A23 - A32) (24) 

V^ = \ (As - A3 + A23 - A32) (25) 

V H = ^= (2e - 2 A2 + As + As - A23 - A32) (26) 



Here the indexes ip and (j) correspond to the notations for radial functions used in Eq.(T12j). 
Each of these projectors produces matrix elements between the two radial basis functions 
from Eq. (|12p with the correct permutation symmetry. Note that the two projectors V^ 
and V^ written above coincide with each other. It can also be shown that the three 
projectors V^,V^ and V^ are orthogonal to each other. In actual computations only the 
upper triangles of the Hamiltonian and overlap matrices are used. Therefore, only the three 
projectors V^,Vf^ and V^ are important in computations of the bound doublet 2 S— states 
in all three-electron atomic systems. 

The approach described above allows one to construct the spatial parts of the total 
variational wave functions with the correct permutation symmetry between all identical 
particles in the three-electron atomic system. In our previous work we have also found 
(see, e.g., jsj) that the same approach works perfectly for all four-, five- and many- 
electron systems. Moreover, the symmetry of the electron spin functions can also be different, 
e.g., for the singlet and triplet states in four-electron systems. Below, the variational wave 
functions for the singlet and triplet states in the four-electron atomic systems are explicitly 
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constructed. The explicit formulas for the spatial parts of trial wave functions are derived 
with the use of the corresponding spatial projectors. 

B. Four-electron atomic systems. Singlet states. 

Numerical computations of the bound states in four-electron atomic systems include the 
non-trivial step of antisymmetrization of all electronic variables, i.e. variables 1, 2, 3 and 4 
in the trial wave function The variational wave function \1/ of any sing let 1 S(L = 0)-state 
in four-electron atomic systems is represented in the form of Eq. fl27|) . i.e. 

^l=o = ipL=o(A; {r i:j })(a(3a(3 + (3a(3a - (3aaf3 - a(3(3a) (27) 
+c/) L=0 (B; {rij})(2aaf3(3 + 2f3[3aa — f3aa(3 — a(3/3a — f3af3a — aj3aj3) 

where if)L =0 (A; {r^}) and (pL=o(B; {r^}) are the two independent spatial parts (= radial 
parts) of the total four-electron wave function. The notations A and B mean that the two 
sets of non-linear parameters associated with if) and <fi are optimized independently. Such a 
trial wave function must be antisymmetric upon all electron variables (or variables 1, 2, 3 
and 4 in our notations), i.e. A e ^ = where 

A e = e — P\2 — Prs — P23 — Pu — -P24 ~ -P34 + -P123 + -P132 + -P124 + -P142 + -P134 (28) 
+-P143 + -P234 + -P243 — -P1234 — -P1243 — -P1324 — -P1342 — -P1423 — -P1432 + P12 ' -P34 

+-P13 ' -P24 + Pu ' -P23 

is the complete four-particle antisymmetrizer. Here e is the identity permutation, while Pij 
is the permutation of the particles i and j. Analogously, the operators P^ and Pijki are the 
permutations of the particles i, j, k and i, j, k, I, respectively. 

By using the procedure described in the previous Section we can find the explicit formulas 
for the corresponding spatial projectors. In fact, by using the explicit form of the trial wave 
function ^ constructed for singlet states in four-electron systems, Eq. (l27|) . and by integrating 
over electron spin coordinates one finds the four following spatial projectors: 

Vw = ^=(2e + 2A2 - A3 - As - Pu ~ A 4 + 2A4 + 2A2A4 + 2A3A4 + 2P M P 2 {29) 

— Pl23 ~ -Pl32 — -Pl24 — "l42 ~ -Pl34 — -Pl43 — -^234 — -^243 — -Pl234 — -Pl243 

+2A324 — A342 — A432 + 2A423) 
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"P^<j> — J (A3 ~~ A3 ~~ P\4 + P24 + A23 — A32 — A24 + A42 + A34 — A43 — A34 + A43(30) 

+ A234 — A243 — A342 + A432) 

P<M = 47! ( 2S ~ 2pl2 + Pl3 + p23 + PlA + p24 ~ 2pM + 2pl2p ' M + 2pl3p24 + 2 ^ 14 ^ 32 ) 

— -Pl23 — -Pl32 — -Pl24 — -Pl42 — -Pl34 — -Pl43 — -^234 — -P243 + -Pl234 + -Pl243 

— 2 A324 + A342 + A432 — 2A423) 

In reality, since = V^, one needs to use only three such operators V^,V^ and V^. 
The use of these three projectors for matrix elements allows one to produce the matrix 
elements which have the correct permutation structure between all four identical particles 
(electrons). Note that all such matrix elements are computed only between the corresponding 
spatial basis functions and do not include any spin function. The explicit formulas for the 
complete set of singlet spatial projectors for four-electron atomic systems, Eqs. (l2"9]) - (I3"2"j) . 
have not been presented in previous publications. First bound state computations of four- 
electron atomic systems with the use of completely correlated wave functions were performed 
by Sims and Hagstrom [8j. Since then many authors have conducted such calculations for 
singlet states in various four-electron systems (see, e.g., jsj and references therein). 



C. Four-electron atomic systems. Triplet states. 

The trial wave function of the triplet 3 S'-state in the four-electron atomic system can also 
be represented in the form with the two independent spin functions xi — ct(3aa — f3aaa and 
X2 = 2aaf3a — f3aaa — a(3aa. The variational expansion takes the form 

\I/ = -0i =o (A; {rij}) ■ (a(3aa — (3aaa) + 4>l =q (B; {ry}) • (2aa(3a — j3aaa — a(3aa) (33) 

where ipL=o(A; {ry}) and 4>l =q (B; {rij}) are the radial parts (also called the spatial part) 
of the total wave function. Here the notation {r^} designates the complete set of fifteen 
interparticle (spatial) coordinates, while the symbol A and B mean the corresponding sets 
of non-linear parameters. Optimization of non-linear parameters in the A— and/or I?-sets 
is performed independently from each other. The trial wave function, Eq. (l34|) . contains 
two electron spin functions \i an d X2 which correspond to the S = 1 and S z — 1 values, 
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where S and S z are the eigenvalues of the total electron spin and its z— projection, i.e. 
S 2 Xi = S(S + l)xi and S z Xi = S z \i- 

For triplet states by using the explicit form of the xi an d X2 functions one can easily find 
the four other spin functions which correspond to the S = 1 and S z = and S — 1 and 
S z = — 1 values. For instance, in the case of Xi — Xi s Pi n function the two spin functions 
Xi" 1 = af3af3 + a/3f3a — f3aa/3 — f3af3a and Xi~^ = ce(3{3(3 — (3a(3(3 correspond to the S = 
1, S z — and S — 1, S z — — 1 values, respectively. The three spin functions Xi? Xi > Xi ^ 
form a regular triplet of spin functions. An analogous triplet of spin functions can be 
constructed for the X2 spin function. To describe experimental situations with no external 
magnetic field present we need to use the spin functions associated with one value of S z , 
say, Sg = 1. Thus in this work we shall always choose Xi — ( a fi — (3a)aa = a(3aa — (3aaa 
and X2 = 2aa(3a — (3aaa — a(3aa (see Eq. fj33|) above). 

In calculations with such a trial wave function, Eq. (I33j) . one needs to know the explicit 
formulas for all three radial (or spatial) projectors V^, V^{= V^) and V^. However, 
numerical calculations using a wave function with two spin functions are computationally 
intensive and were not attempted in the current study. Accordingly, we did not attempt to 
derive the associated projectors. Instead, we performed some computations of the triplet 
states in four-electron atomic systems with the use of one spin function xi = = 
a(3aa — (3aaa only. The variational expansion Eq.f l33|) is now written in the form 

\& = vpL=o(A; {rij}) ■ (a/3 — (3a)aa = i(jl =q (A] {ry}) ■ (a(3aa — (3aaa) (34) 

Now, we need to obtain the spatial part of the total wave function with the correct permu- 
tation symmetry between all identical particles 1, 2, 3 and 4 (electrons). The corresponding 
spatial projector is obtained in this work by calculating the explicit expression for the fol- 
lowing spin expectation value 

= CY^ s ahcd ((a(3 - (3a)aa \ P abcd \ (a/3 - (3a)aa)P abcd (35) 

abed 

where C is the normalization factor, while the integers s a bcd are defined from the explicit 
form of the complete four-particle (or four-electron) antisymmetrizer A e given in Eq.( l28|) . 
After some algebra one finds the explicit formula for the corresponding spatial projector 

= -^-j={2e + 2A2 — A3 — P23 — Pu ~ P24 ~ 2F34 — 2A2A4 — A23 — A24 (36) 
— A32 — A42 + A34 + A34 + A43 + A43 + A432 + A234 + A243 + A342) 
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This projector creates the spatial part of an arbitrary matrix element needed in bound state 
computations of the triplet 3 S'— states in an arbitrary four-electron atomic system. Such 
a matrix element has the correct permutation symmetry among all four identical particles 
(electrons). Explicit formulas for the spatial projectors which correspond to the triplet states 
have not been published previously. 

By using the formulas presented above one can perform accurate computations of the 
triplet bound states in various four-electron atomic systems. As follows from our results of 
such calculations (see, e.g., the method described above allows one to determine various 
expectation values in four-electron atomic systems to relatively high numerical accuracy. In 
particular, such expectation values can be computed for all positive and negative powers of 
relative coordinates r^. In general, the expectation value of any regular function of the ten 
relative coordinates can be computed to very good numerical accuracy. Analogous ex- 
pectation values which contain delta-functions of the relative coordinates and their products 
with the regular functions of the relative coordinates also do not present any problem for 
numerical computations. Real problems arise in computations of expectation values which 
include products of delta-functions with the corresponding electron spin functions, e.g., 
(aiSm) and {oci^jSm^Nj) > where iV designates the nucleus, while i means i-th electron. The 
first expectation value (cuSm) represents the single-electron density of a-electrons on the 
atomic nucleus. Analogous expectation values can be computed in the case of (3— electrons. 
A very poor convergence of similar expectation values means that another spin function 
must be included in computations. In Eq.(l33l) such a spin function is called the second 
(electron) spin function \2- 

V. GENERALIZATION TO THE FIVE- AND SIX-ELECTRON ATOMIC SYS- 
TEMS. 

The method described above allows one to construct the properly antisymmetrized trial 
functions for three- and four-electron atomic systems. Formally, our method can be gener- 
alized to the five, six and many-electron atomic systems. However, its direct generalization 
is very difficult, since the proper antysymmetrization of the trial basis functions and their 
linear combinations becomes extremely difficult in the case of many-electron atoms with 
A — 1 > 5 electrons, i.e. for A > 6. Here and everywhere below, A is the total number of 
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particles in the system, while A — 1 is the total number of electrons. This means to construct 
the trial wave functions for five-, six- and many-electron atoms one needs to apply A\(> 120) 
different interparticle permutations to the non-symmetrized basis function. The presence 
of a very large number of terms in each wave function drastically complicates the explicit 
expressions for the spatial projectors mentioned above. For instance, such a spatial projec- 
tor constructed for the B-atom (five-electron atom) must include 120 different terms. Some 
of these terms can be equal zero identically, but in any case the total number of remaining 
terms is still very large. Therefore, it is important to develop some effective methods which 
can be used to operate with a very large number of terms in the trial wave functions. Our 
current hopes rely on the two following methods. The first approach is based on the use of 



16| . In this approach all 



various symbolic-algebra computational platforms such as Maple 
integrations over spin variables can be performed analytically. The expressions for all spatial 
projectors are never written explicitly, but they are used internally by this computational 
platform. Note also that for some basis sets the action of any interparticle permutation P a bc... 
on the basis wave functions is reduced to the permutation of the corresponding non-linear 
parameters in these functions. In particular, this is the case for variational expansion defined 
by Eqs.(jl]) and ([5]). This means that actual permutations of the non-linear parameters in 
the basis wave functions, Eqs.Q and ([5]), can always be applied instead of the permutation 
of the relative coordinates. This drastically simplifies the explicit construction of the com- 
pletely symmetrized trial wave functions. The permutation of the non-linear parameters in 
the basis wave functions can be combined with the analytical integration over spin variables 
in the total wave function. This can be used in the future methods. 

The second approach is based on the relations which exist for the spin functions in 
three-, four- and many-electron systems. For instance, the second spin function \2 used in 
numerical computations of the triplet states of four-electron systems (see above) is obtained 
from the X2 spin function known for the doublet states in three-electron systems. Formally, 
we can write X2(l, 2, 3, 4) = %2(1, 2, 3)a(4), where a(4) is the spin function of the additional 
(= fourth) electron, while symbols (1,2,3,4) and (1,2,3) designate the systems of four and 
three electrons, respectively. The notation x 2 (l>2,3) stands for the second spin function of 
the doublet state in three-electron system, while X2(l>2,3,4) means the second triplet spin 
function of the four-electron system. Analogous relation exists between another (first) triplet 
spin function of the four-electron atomic systems, xi(l,2,3,4), and xi(l,2,3) used above 
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(see, e.g., Eq.( TT2]) ) for the doublet states in three-electron atoms/ions. By studying this and 
other similar relations between spin functions we can find some useful connections between 
the spatial projectors constructed for three- and four-electron systems. This approach can 
also simplify methods and algorithms which must be developed in the future for systems 5 
or more electrons. 

VI. NUMERICAL RESULTS 

To illustrate our method in applications to actual three- and four-electron atomic systems 
let us briefly describe the results of variational computations of bound states in the three- 
electron Be + ion in its l 2 S'-state and four-electron Be atom in its l 1 ^— and 2 3 S'-states. For 
simplicity all nuclear masses were assumed to be infinite in such calculations. A separate 
group of calculations have been performed for the 2 3 5'— electron state in the six-body oxygen- 
muonic ion 8+ fi~e^ , where /i~ is the negatively charged muon. This positively charged ion 
(q = + 3) is a well bound atomic system which contains the composite 'nucleus' (0 8+ + //~) 
with overall 'nuclear' charge +7 and four atomic electrons. Below, we consider the 16 
nucleus only. In our calculations of the 8+ u~e7 ion we used M = 29156.9457 m e for the 

n n 

nuclear mass of the oxygen- 16 nucleus and = 206.768262 m e [14j, UM- 

Numerical results of our computations can be found in Table I where one finds the total 
energies E and some other bound state properties expressed in atomic units, where h = 1, m e 
= 1 and e = 1. The electron state of each atomic system is shown in the following brackets. 
In the case of the 8+ ^x~eJ ion the notation 2 3 S' e stands for the 2 3 S , e -triplet electron state 
in this system. The muonic quasi-nucleus 8+ + /i~ is in its ground l 1 ^— state. This is 
always assumed, but not shown in our notation. Note also that for bound state properties 
presented in Table I the index e means the electron, while the index N denotes the atomic 
nucleus. For each energy shown in Table I only 9 decimal digits are presented. In general, 
optimization of the non-linear parameters in variational expansion Eqs.()l]) and fl5]) always 
decreases the total energies. On the other hand, small variations in a few last decimal digits 
are not critically important for our present purposes. 

As follows from Table I our method provides very good numerical accuracy for doublet 
states in three-electron atoms and ions. On the other hand, this method also works perfectly 
for singlet and triplet four-electron atomic systems. It is very likely that the analogous 
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procedure can be developed for five-, six- and many-electron atomic systems. However, for 
atomic systems with five and more electrons one finds a number of additional problems and 
direct generalization of our method is very difficult (see discussion in the previous Section). 
Note also that the overall convergence rates of the radial variational expansion, Eqs.(|3J) - 
for three- and four-electron atomic systems are comparable with each other. It seems to 
be very strange, but we need to remember that the number of non-linear parameters in each 
basis function rapidly increases as the total number of bodies A in system increases. For 
four-electron atomic systems each of these basis function contains 10 non-linear parameters 
(each of them varied independently), while for three-electron atomic systems one finds only 6 
such parameters in each basis function. Briefly, this means that the overall 'flexibility' of the 
four-electron trial function is comparable with the analogous 'flexibility' of the three-electron 
trial functions. 

VII. CONCLUSION 

We have considered the problem of accurate computations of bound states in three- and 
four-electron atomic systems. The method developed in this study allows one to construct 
variational wave functions for an arbitrary bound state in three- and four-electron atomic 
systems. All such trial wave functions have the correct permutation symmetry (with respect 
to all permutations of identical particles, i.e. electrons). It is important to note that in our 
method the total number of independent spin functions can be varied. Numerical compu- 
tations can start with the use of one electron spin function only. The second spin function 
can be introduced later to improve the overall convergence of the results. Our procedure 
can be generalized to bound states with the non-zero angular momentum L. Such a gener- 
alization is straightforward, but it requires extensive use of additional notations, application 
of special methods developed in theory of angular momentum and substantial explanations. 
Variational calculations of the bound P(L = 1)— states in five-electron atomic systems will 
be considered in our next study. The ground state of the B-atom is the P(L = 1)— state. 

Our method is based on explicit constructions of the total wave functions for various 
bound states in three- and four-electron atomic systems. The unified procedure has been 
applied to an example of each type of system. The central part of the procedure is the 
construction of spatial projectors with the correct permutation symmetry between all iden- 
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tical particles (electrons). This method was originally developed for three-electron atomic 
systems by Larsson [9J (see also [lOj). We generalized this procedure to the cases when 
a number of different spin functions are used in computations. In addition, we have con- 
structed spatial projectors needed in calculations of the singlet and triplet bound states in 
four-electron atomic systems. Currently, the total energies and other properties of bound 
states in four-electron atomic systems can be determined to the accuracy which is better 
han the accuracy of old bound state calculations performed for two-electron atomic systems 



171 ]. In general, by using various optimization strategies for non-linear parameters from the 



trial wave functions one can obtain very accurate variational energies and highly accurate 
wave functions. Such wave functions can be used in the following accurate computations of 
different bound state properties, including various relativistic and QED corrections. 



[1] L.D. Landau and E.M. Lifshitz, Quantum Mechanics. Non- Relativistic Theory (3rd. ed., Ox- 
ford, England, Pergamon Press (1977)). 
[2] A.M. Frolov and D.M. Wardlaw, Phys. Rev. A 79, 064501 (2009). 
[3] F.G. Galvez and E. Buendia, J. Chem. Phys. 117, 6071 (2002). 
[4] C. Froese Fisher, J. Phys. B 26, 855 (1993). 

[5] A.M. Frolov and D.M. Wardlaw, Phys. Rev. A 78, 042506 (2008). 

[6] N.N. Kolesnikov and V.I. Tarasov, Yad. Fiz. 35, 609 (1982), [Sov. J. Nucl. Phys. 35, 354 
(1982)]. 

[7] A.M. Frolov and D.M. Wardlaw, JETP 135, 583 (2009). 
[8] J.S. Sims and S.A. Hagstrom, Phys. Rev. A 4, 908 (1971). 
[9] S. Larsson, Phys. Rev. 169, 49 (1968). 
[10] FW. King, Phys. Rev. A 40, 1735 (1989). 

[11] M.I. Petrachen and E.D. Trifonov, Applications of Group Theory in Quantum Mechanics 

(Dover Publications, Inc., New York, 2009), Chpts. XV - XIX. 
[12] V. Heine, Group Theory in Quantum Mechanics (Dover Publications, Inc., New York, 1993), 

Chpts. II - IV. 

[13] P.R. Halmos, Finite- Dimensional Vector Spaces, (2nd. ed., D. Van Nostrand Company Inc., 
New York, (1958)), Clip. III. 



20 



[14] The NIST Reference on Constants, Units and Uncertainty, see: 

http:/ / physics.nist.gov/cuu/Constants/index.html 
[15] CRC Handbook of Chemistry and Physics, 85th Edition, Ed. D.R. Lide, (CRC Press, Inc., 

Boca Raton, Florida, 2004). 
[16] A product of Waterloo Maple Inc., Waterloo, Ontario, Canada (see: 

|http://www.maplesoft.com[ ). 
[17] H.A. Bethe and E.E. Salpeter, Quantum Mechanics of One- and Two-Electron Atoms 

(Springer- Ver lag, Berlin, 1957). 



21 



TABLE I: The non-relativistic energies and other properties determined for the S— states in some 
three- and four-electron atoms, ions and muonic ions (in atomic units). 





Be+ (1 2 S) 


Be (l 1 ^) 


Be (2 3 S) 


8+ /x-eJ (2 3 S e ) 


E 


-14.3247627 


-14.6673323 


-14.4300595 


-6619.33457 


Kn) 


2.65796 


2.10684 


2.03603 


4.77655 


(Tee 1 ) 


1.08200 


0.72912 


0.61933 


1.26072 


{r e N) 


1.03379 


1.49297 


2.63085 


0.98798 


(r ee ) 


1.75565 


2.54516 


4.70847 


2.43483 
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